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Abstract—In this paper, a general lumped parameter bear- 
ing current model is presented. The lumped parameters are 
derived from the analytical descriptions of the bearing current 
paths and the numerical field simulations based on the finite 
element method. The model considers frequency dependence. 
The resulting model delivers the common mode current and 
the circulating bearing current is determined by modeling the 
magnetic coupling between the common mode current path and 
the circulating bearing current path. The derived circuit can 
be applied to AC machines and is verified by measurements 
of a permanent magnet excited synchronous machine, which is 
supplied by an off-the-shelf motor drive inverter. 


Index Terms—Bearing currents, circulating bearing current, 
common mode current, common mode voltage, magnetic cou- 
pling, transmission ratio 


I. INTRODUCTION 


Although, the bearing current issue in electrical machines 
is known for nearly 100 years [1], it is still a relevant ongoing 
research topic. Failures caused by bearing currents lead to 
heavy mechanical damages in electrical machines resulting 
in high maintenance costs [2]. 

Pulse width modulated (PWM) drives, which use very fast 
switching devices such as silicon carbide (SiC) or gallium 
nitride (GaN) switches, in order to reduce the converter losses, 
may even worsen the issue [3, 4]. These fast switching devices 
generate a common mode (CM) voltage with a very high 
du/dt, which leads to a CM-current with very high fre- 
quency (HF) components [5]. The HF-CM-current is mainly 
conducted through the stator-winding stator-core capacitance, 
back to the source through the grounded machine housing (see 
path | in Fig. 1). This causes the so-called circulating bearing 
(CB) current [6-8]. The CB-current is an HF inductively 
coupled current, which arises in the conductive loop formed 
by the stator core-housing-non drive end (NDE) bearing-shaft- 
rotor core-drive end (DE) bearing-housing (see path 2 in 
Fig. 1). It is worth mentioning that there also exists other types 
of bearing currents, such as the so-called electric discharge 
machining (EDM) bearing currents or du/dt-currents (see 
path 3 in Fig. 1). The EDM-currents are capacitively coupled 
from the stator winding to the bearings. They are conducted 
when the voltage across one of the bearings exceeds a certain 
limit and a breakdown occurs. The du/dt currents are also 


978-1-5386-6376-9/19/$3 1.00 ©2019 IEEE 


Djordje Vukovic, Sebastian Blume 
B&R Industrie-Automation AG 
CH-8500 Frauenfeld, Switzerland 


Housing (H) 


g < 
Pathl lem 


(2) 
€ 


Stator end 
winding 


SEW-H Stator winding 
star point 


¥. lop 
[Path 2 


Isolation 
layer 


Bearing (DE) 
Rotor 
core (RC) 


Symmetry axis 


Fig. 1. Simplified representation of the CM-current (icy, icm s path 1), 
the CB-current (icp, path 2) and the EDM-currents (tgpm, tdu/at, paths 3) 
paths and in an AC-machine. The excitation source of these currents is the 
CM-voltage wom. The definitions of the abbreviations are listed in Tab. I. 


capacitively coupled from the stator winding to the bearings. 
They are flowing as a displacement current through the 
bearing capacitance [9, 10]. However, these currents are not 
discussed in this paper. 

In the current literature, the bearing currents are determined 
with the help of equivalent lumped RLC-circuit models, 
which represent the equivalent bearing current paths in the 
electrical machine [5—13]. The RLC-parameters are derived 
from analytical approximations of the machine geometry or 
finite element method (FEM) based simulations [14, 15]. 
However, the magnetic interaction between the CM-current 
and the CB-current is highly simplified in the existing models 
[6, 8]. Therefore, in this paper, an improved lumped parameter 
model, which accurately considers the magnetic coupling 
between the CM-current and the CB-current is presented 
and verified by measurements of a permanent magnetic syn- 
chronous machine (PMSM). The calculations of all frequency 
dependent lumped circuit elements, as well as the FEM based 
derivation of the improved magnetic coupling parameters are 
given. 

This paper is organized as follows. In section II, the general 
bearing current model is introduced and all governing equa- 
tions for calculating the CM- and CB-currents are derived. 
The derivation of the lumped RLC-parameters is given in 
section III, where the main part of this section is concentrating 
on the determination of the magnetic coupling parameters. 
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TABLE I 
DESCRIPTION OF THE USED ABBREVIATIONS IN FIG. 1 AND FIG. 2 


Abbreviation Description 

SW Stator winding 

EW End winding region 

SWA Part of the stator winding inside the stator core 
SEW Part of the stator winding outside the stator core 
SC Stator core 

RC Rotor core 

H Housing 

B(NDE) Bearing non drive end 

B(DE) Bearing drive end 

SFT Shaft 

Cab Grounding cable 


Finally, in section IV the measured CM- and CB-currents are 
verified with the currents delivered from the proposed model. 


II. GENERAL COMMON MODE AND BEARING CURRENT 
MODEL 


Figure 2 (a) shows the equivalent lumped parameter model 
according to Fig. 1. There, path 1 is related to the CM-currents 
icm and icy’, path 2 is related to the CB-current icg and 
path 3 is related to the EDM-current igpy and du/dt-current 
dawat- Since, the capacitor Cswa-sc is significantly bigger 
than the other capacitors in the circuit (compare capacitances 
in Tab. II), it can be safely assumed icy œ% icm’. The 
bearing resistance Rpnpe)pE), aS well as the non linear 
impedance Zgane) oe) represent the combined values of the 
bearing from the NDE-side and the DE-side. By assuming a 
breakdown of both bearings at the same time (Agnpe),(DE) © 
0, Zenne), oE) © 0), the CB-current becomes the dominating 
current through the bearings [9] and the model can be reduced 
to the model depicted in Fig. 2 (b). In this case, also path 3 
is neglected. An even more compact model is shown in 
Fig. 2 (c). Due to the frequency dependence of the lumped 
circuit elements, the following analysis is carried out in 
the complex frequency domain to determine the CB-current. 
The currents Jom(jw) and Icg(jw), as well as the voltages 
Ucm(jw) and Ucg(jw) are the complex representations of 
their time domain signals icm, icp and ucm, Ucp, respectively 
w is the angular frequency. After solving the model in the 
frequency domain, the spectral components of the currents 
are finally transformed back into the time domain by using 
the well known fourier series [16]. 

Applying Ohm’s law on the reduced bearing current model 
in Fig. 2 (b) leads to 


; Ucm (jw) 
I = 1 
cm(jw) Zai (1) 
where Zeaq iS given as 
Zred (Jw) =Zeq (jw) +Rsw(w)+A(jw)+Ru(w)+ 
(2) 
+ Reav(w) F jwLcar(w) 
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Fig. 2. (a) Equivalent lumped parameter model according to Fig. 1. Path 1 
is related to the CM-current icp, path 3 is related to the circulating bearing 
current icp and path 3 is related to the EDM- and du/dt-currents igpm 
and iguyat- (b) Reduced equivalent lumped parameter model and (c) compact 
controlled source lumped parameter model, which is equivalent to (b) and 
used for the verification in section IV. 


and 


1 
jw(Cswa-sc + Csew-sc + Csew-n) 


A(jw) = (3) 


The governing equations for the magnetic coupling between 
path 1 and path 2 are obtained as 


Ucsi (jw) = jwLsw(w)om (jw) —jwM (w) Ice (jw) (4) 
Uc2(jw) = jwM (w)Icm (jw) — jwLe(w) Ics (jw) (5) 


where, M (w) is the frequency dependent mutual inductance, 
Lsw(w) is the frequency dependent self inductance of the 
stator winding and L2(w) is the frequency dependent self 


inductance of the CB-current path 2, which is given by 
Lo (w) = Lsc(w) + Ire (w) + Lew (w) . (6) 
By substituting [cg with 
Uce2(jw) 
I SA E T 
cB (jw) Ra (7) 
226 
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TABLE I 
EQUIVALENT CIRCUIT CAPACITANCES 


Location Capacitance (pF) FEM method 
CsEw-sc 209.71 3D 
CsEW-RC 6.00 3D 
Csew-H 113.04 3D 
CsEW-B(NDE),(DE) 3.71 3D 
CsEW-SFT 11.81 3D 
Cywa-sc 46 407.77 2D 
Cswa-RC 3.50 2D 
in (4) and (5), leads to 
Í : . : Ui ju 
Ue (jus) =e Lsw (wo) Teo (jer) -juM E 
Rce (w) 
. jwM(w)Rcg(w) ; 
U, = I 9 
cB2(jw) ney alae) cm (jw) (9) 
where 
Res (w) = Rsc(w) + Rrc(w) + Rew(w) (10) 
Inserting (9) in (8), results in 
Ucai (jw) = Tem (jw) Zeg (Jw) (11) 
where 
w? M? (w) 
Zealjw) = jwL H - . 12 


Finally, by the substituting (9) in (7), the complex transfer 
characteristic of cg is obtained as 


Icg(jw) = juM(w) 
Iem(jw) = Reg (w) + jwLe(w) 


Gigs (jw) = (13) 
with 

M(w) = wepcaveL2(w) (14) 
and wepcavg is the magnetic transmission ratio. Notice: Equa- 
tions (1) and (13) are represented by the model in Fig. 2 (c). 
The derivations of the lumped frequency dependent parame- 
ters of the reduced bearing current model is presented below. 


III. DERIVATION OF THE CM- AND CB-CURRENT 
CIRCUIT LUMPED PARAMETERS 


A. Derivation of the RLC-circuit elements 


The FEM based derivation of the capacitances, as well as 
the resistance and the inductance of the end winding region 
Lew and Rew, is performed in the same way as presented 
in detail in [8]. The resistive and inductive stator core and 
rotor core parameters Rsc, Rec, Lsc and Lec are analytically 
obtained according to [17]. The derived values are listed in 
Tab. II, respectively shown in Fig. 3. In addition, the measured 
CM-impedance Zcm meas is depicted in Fig. 3. 
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Fig. 3. Finite element method obtained (Lew, Rew), analytically ob- 


tained (Lsc, Lac, Rsc, Rec) and measured (ZcM meas) frequency dependent 
lumped parameters of the reduced bearing current model. 


B. Derivation of the magnetic coupling parameters 


To obtain the mutual inductance M/(w) (see (14)), only 
an eddy-current analysis would be not sufficient, since the 
source current of this model flows as a displacement current 
through a capacitive path (the capacitive coupling between 
the stator winding and the stator core). Thus, a Full-Maxwell 
field analysis including inductive and capacitive couplings 
should be performed. This analysis would be too demanding 
in terms of CPU time and memory and has considerable 
stability problems as the considered frequency range goes 
from RF- (MHz) towards low frequencies (Hz). To avoid the 
mentioned difficulties the following method is applied in this 
work. 

Starting from the second Maxwell equation in its complex 
form [18] 


— = 2 
Vx H=J+jwD (15) 
and taking the divergence on both sides 
— d _2 
V-(VxH)=V-J+V-jwD (16) 
results in 5 - 
0O=V-J+V-jwD (17) 
and finally in 
> > 


where H is the magnetic field and D is the electrical flux 
density. After introducing the material relations, (18) can be 
written as 
=> 2 

V-Eo=—V-jwE e (19) 
where E is the electrical field, o is the electrical conductivity 
and e is the electrical permittivity. By comparing the material 
coefficients results in 


lolpc ==> |w- €lac (20) 
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Fig. 4. A 10° section of the PMSM is modeled in 3D. (a) View 1 of the 3D-SCD-FEM model of the investigated PMSM. The stator coil is shown within 
a breakout view of the stator core. (b) View 2 of the 3D-SCD-FEM model of the investigated PMSM. Only the Stator coil (orange), the stator core (grey), 
the insulation layer (yellow) and a part of the housing (green) is needed for solving the SCD-BVP. (c) Axis symmetric 3D-FEM model, where all parts of 


the PMSM are modeled for solving the MS-BVP. 


By using the interpretation (20) the capacitive displacement 
current can be replaced with equivalent conduction current 
that has an identical path and the same current density. Having 
this mathematical conversion a 10° section of the PMSM is 
modeled in 3D. The current density distribution of this current 
source can be determined by solving the three dimensional 
stationary current distribution (SCD) boundary value problem 
(BVP). The 3D-SCD-BVP is defined in the computational 
domain Q; as [19] 


V- (opc: Vy) =0 in Qı CR? (21) 
Y = 0 on Op1 94 (22) 
n opc: Vg=0 on AnD, (23) 

I 
n -Vy = —<— on M (24) 

ODC S 
and finally, solving for the current density results in 
> 

Jsw=—opcVy with opc =w : eac. (25) 


The declarations of the 3D-SCD-BVP are depicted in 
Fig. 4 (a) and (b). The DC-current Ipc is applied at the surface 
S of the inlet (boundary n291) and the scalar electrical 
potential ~ is applied at the top of the outlet (boundary 
0p, 1). The boundary Oy;2, (blues lines) is applied on all 
conducting outer surfaces, which means, no current density 
component perpendicular to the conducting surfaces exists. 
In order to solve the 3D-SCD-BVP in the Qı domain, only 
the stator coil, the stator core and parts of the housing are 
modeled, as can be seen in Fig. 4 (a). 

The mutual inductance Mpc, as well as the CB-current loop 
inductance Lə pc are obtained by solving the 3D magneto- 
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static (MS)-BVP, which is given in the computational domain 
Q as [19] 


il > > 
Vx [ev Al =Jsw in QCR? (26) 
H 
y > 
n x A=0 on OpQ (27) 
ae > 
n xVx A=0 on ON? (28) 


where A is the magnetic vector potential and u is the magnetic 
permeability. At the axis of symmetry, only a tangential 
component of the vector potential is non-zero (boundary 
OpQ) and there exists no magnetic field outside the air box 
(boundary OyQ), as shown in Fig. 4 (c). For solving the 3D 
MS-BVP in the Q domain, all parts of the PMSM are included 
(see Fig. 4 (c)). Finally, the magnetic transmission ratio wepc 
is calculated as [20] 


Mpc(w) 


~ Depew)’ 
The close invariance in frequency of the magnetic coupling 
parameters is clearly shown in Fig. 5. Thus, the mean value 
of uepc(w) (see Tab. IM) is used for calculating M (w) in (14) 
with 


wepc(w) (29) 


UEDC, avg = Mean(uepc(w)) (30) 


IV. MEASUREMENT AND SIMULATION RESULTS 


In the following, the simulated CM- and CB- 
currents are verified with measurements on a PMSM 
(8LSP96.R2015P644-3, [21]), which is driven by a inverter 
based motor drive (ACS880-01-061A-3, [22]). The main 
parameters of the motor drive system are listed in Tab. IV 
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TABLE III 
AVERAGED MAGNETIC COUPLING PARAMETERS 


Parameter Value 
L2 peave 2.36 x 1071 H 
Mpcave 3.85 x 1072 H 
UEDC,avg 1.64 x 10-1 
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Fig. 5. FEM obtained mutual inductance Mpc, self inductance of the CB- 


current loop Lo pc (see Fig. 4 (c)) and the magnetic transmission ratio wepc. 


and the measurement circuit is depicted in Fig. 6. The 
investigated PMSM is equipped with a hybrid bearing on 
its NDE-side, which is bypassed with a copper stripe for 
measuring the CB-current iCB,meas. The CB-current probe 
(CPO31A), as well as the CM-current probe (CP031), [23] 
have a bandwidth of 100 MHz respectively a sensitivity of 
1mA/div and 10mA/div. Further, the test machine has 
been mounted isolated from the ground and only a defined 
single point at the housing provides the grounding. The shaft 
voltage Usrrmeas is Measured on the DE-side with respect 
to the housing and the measured CM-voltage UCM,meas iS 
obtained by measuring all phase-to-ground voltages and 
calculating the mean value. Further, the phase-to-ground 
voltages are measured as close as possible to the machine in 
order to exclude the influence of the supply- and grounding 
cables on the measurement of UCM,meas. The machine is 
operated under no-load conditions. 

For verifying the proposed CB-model, the circuit in Fig. 2 (c) 
is used and Zeq is replaced by the measured CM-mode 
impedance Zcm meas (See Fig. 3), which is obtained by 
electrically connecting the endings of the stator winding and 
measuring the impedance related to the grounded point of 
the housing. The depicted period of ucmmeas in Fig. 7 (a) is 
transformed into its spectral components and input into the 
frequency depended model. The considered frequency range 
is between 30 Hz and 5 MHz. 

Figure 7 shows all measured values according to Fig. 6 and 
the simulated CM-current zcmsim, as well as the simulated 
CB-current icgsim. Figure 8 shows the waveforms of the 
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TABLE IV 
DATA SHEET PARAMETERS OF THE INVESTIGATED PMSM [21] AND THE 
USED MOTOR DRIVE INVERTER [22] 


Parameters PMSM 


Nominal speed ny | (rpm) 1500 
Nominal torque My | (Nm) 350 
Nominal power Py (W) 54978 
Nominal current (RMS) In (A) 107.3 
Parameters motor drive inverter 
Nominal power Py (W) 30 000 
Nominal current (RMS) Ix (A) 61 
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Fig. 6. Measurement setup for measuring the CM-current iCM,meas, the CB- 
current 2cB meas. the CM-voltage UCM,meas and the shaft voltage usrrmeas- 


zoom area in Fig. 7. The time range is exemplary chosen 
where the shaft voltage is low and a maximal CB-current 
occurs. By comparing the measured and the simulated first 
maximal peak values of the CM-currents in Fig. 8 (c) (see 
regions Yı and Yə) the relative error is 14.7% for the 
negative peak and 11.7% for the positive peak. The relative 
error in case of the CB-currents in Fig. 8 (d) (see regions 
X, and X2) is 38.6% for the negative peak and 15.8% for 
the positive peak. The higher deviation of the CB-current 
can be explained by the fact that the bearing resistance is 
assumed to be zero during the breakdown and the impedance 
of the additional copper stripe in the bearing current circuit 
is neglected as well. Therefore, the assumed impedance 
of the CB-current circuit is too low and the calculated 
CB-current is too high, which results in a higher deviation. 
However, the proposed model shows good agreement with 
the measurement results in time response. 


V. CONCLUSION 


In this paper, an improved lumped parameter circulating 
bearing current model is presented and verified by measure- 
ments on a permanent magnet excited synchronous machine. 
The derivations of the frequency dependent lumped param- 
eters are given and a detailed finite element based method 
for deriving the magnetic transmission coefficient from the 
geometry of the electrical machine is carried out. The relative 
error between the peaks of the circulating bearing current 
of the proposed model and measurements is in the range 
of 15.8% and 38.6%, whereas the waveforms show good 
agreement in time response. 
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The peaks in the regions Y1, Y2, Xı and Xə are used for the evaluation. 


230 


[1] 


[4] 


11] 


23 


REFERENCES 


P. L. Alger and H. W. Samson, “Shaft currents in electric machines,” 
Journal of the American Institute of Electrical Engineers, vol. 42, 
no. 12, pp. 1325-1334, Dec 1923. 

RENOWN-ELECTRIC. (2018) Avoiding Motor Shaft Voltage & 
Bearing Current Damage. [Online]. Available: https://www.magnetec. 
de/fileadmin/pdf/bearing-current-damage- guide.pdf 

D. Han, S. Li, W. Lee, W. Choi, and B. Sarlioglu, “Trade-off between 
switching loss and common mode EMI generation of GaN devices- 
analysis and solution,” in JEEE Appl. Power Electron. Conf. Expo. 
(APEC), March 2017, pp. 843-847. 

K. Euerle, K. Iyer, E. Severson, R. Baranwal, S. Tewari, and N. Mohan, 
“A compact active filter to eliminate common-mode voltage in a SiC- 
based motor drive,’ in JEEE Energy Convers. Congr. Expo. (ECCE), 
Sept 2016, pp. 1-8. 

S. Ogasawara and H. Akagi, “Modeling and damping of high-frequency 
leakage currents in PWM inverter-fed AC motor drive systems,” IEEE 
Trans. Ind. Appl., vol. 32, no. 5, pp. 1105-1114, Sept 1996. 

A. Muetze and A. Binder, “Calculation of Circulating Bearing Currents 
in Machines of Inverter-Based Drive Systems,’ JEEE Trans. Ind. 
Electron., vol. 54, no. 2, pp. 932-938, April 2007. 

S. Chen and T. A. Lipo, “Circulating type motor bearing current in 
inverter drives,’ IEEE Ind. Appl. Magazine, vol. 4, no. 1, pp. 32-38, 
Jan 1998. 

C. Jaeger, I. Grinbaum, and J. Smajic, “Numerical Simulation and 
Measurement of Common-Mode and Circulating Bearing Currents,” in 
22nd. Int. Conf. Elect. Mach. (ICEM), Sept 2016, pp. 486-491. 

D. Schroder, Leistungselektronische Schaltungen: Funktion, Auslegung 
und Anwendung. Springer-Verlag, 2008. 

T. Plazenet, T. Boileau, C. Caironi, and B. Nahid-Mobarakeh, “An 
overview of shaft voltages and bearing currents in rotating machines,” 
in IEEE Ind. Appl. Soc. Annu. Meeting, Oct 2016, pp. 1-8. 

Fei Wang, “Motor Shaft Voltages and Bearing Currents and Their Re- 
duction in Multilevel Medium-Voltage PWM Voltage-Source-Inverter 
Drive Applications,” IEEE Trans. Ind. Appl., vol. 36, no. 5, pp. 1336- 
1341, Sept 2000. 

S. Chen, T. A. Lipo, and D. Fitzgerald, “Modeling of Motor Bearing 
Currents in PWM Inverter Drives,’ IEEE Trans. Ind. Appl., vol. 32, 
no. 6, pp. 1365-1370, Nov 1996. 

D. F. Busse, J. M. Erdman, R. J. Kerkman, D. W. Schlegel, and 
G. L. Skibinski, “The Effects of PWM Voltage Source Inverters on 
the Mechanical Performance of Rolling Bearings,’ IEEE Trans. Ind. 
Appl., vol. 33, no. 2, pp. 567-576, Mar 1997. 


] A. Bubert, J. Zhang, and R. W. D. Doncker, “Modeling and measure- 


ment of capacitive and inductive bearing current in electrical machines,” 
in Brazilian Power Electron. Conf: (COBEP), Nov 2017, pp. 1-6. 

J. hoon Jun, C. kyu Lee, and B. il Kwon, “The analysis of bearing 
current using common mode equivalent circuit parameters by FEM,” in 
Int. Conf. Elect. Mach. Syst., vol. 1, Sept 2005, pp. 49-51 Vol. 1. 

I. N. Bronstein, J. Hromkovic, B. Luderer, H.-R. Schwarz, J. Blath, 
A. Schied, S. Dempe, G. Wanka, and S. Gottwald, Taschenbuch der 
Mathematik. Springer-Verlag, 2012, vol. 1. 

P. Maki-Ontto et al., “Modeling and reduction of shaft voltages in 
AC motors fed by frequency converters,’ Ph.D. dissertation, Helsinki 
University of Technology; Teknillinen korkeakoulu, 2006. 

P. Leuchtmann, Einführung in die elektromagnetische Feldtheorie. 
Pearson Studium Miinchen, 2005. 

J. Smajic and National Agency for Finite Element Methods and 
Standards, How to Perform Electromagnetic Finite Element Analysis. 
NAFEMS, 2016. 

M. Albach, Grundlagen der Elektrotechnik 1. 
Miinchen, 2004. 

B&R—Industrial Automation. (2018) PMSM technical data. 
[Online]. Available: https://www.br-automation.com/de-ch/produkte/ 
antriebstechnik/synchronmotoren- 81s/kuehlart- p/baugroesse-9/ 
8lsp96ee01 5ffgg-3/ 

ABB Ltd. (2018) ACS880_Catalog_3AUA0000118315_RevN_DE(3). 
[Online]. Available: https://new.abb.com/drives/de/ 
frequenzumrichter-fuer-niederspannung/industrial-drives/ 
acs880-single-drive-frequenzumrichter/acs880-01 

Teledyne Lecroy. (2016) Current Probes. [Online]. Available: http: 
//cdn.teledynelecroy.com/files/pdf/current- probes-ds.pdf 


Pearson Studium 


Authorized licensed use limited to: b-on: Universidade de Coimbra. Downloaded on December 20,2023 at 13:12:56 UTC from IEEE Xplore. Restrictions apply. 


